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Abstract. We present a classical linear response theory for a magneto-dielectric 
material and determine the polariton dispersion relations. The electromagnetic field 
fluctuation spectra are obtained and polariton sum rules for their optical parameters 
are presented. The electromagnetic field for systems with multiple polariton branches 
is quantised in 3 dimensions and field operators are converted to 1-dimensional forms 
appropriate for parallel light beams. We show that the field-operator commutation 
relations agree with previous calculations that ignored polariton effects. The Abraham 
(kinetic) and Minkowski (canonical) momentum operators are introduced and their 
corresponding single-photon momenta are identified. The commutation relations 
of these and of their angular analogues support the identification, in particular, of 
the Minkowski momentum with the canonical momentum of the light. We exploit 
the Heaviside-Larmor symmetry of Maxwell’s equations to obtain, very directly, the 
Einsetin-Laub force density for action on a magneto-dielectric. The surface and bulk 
contributions to the radiation pressure are calculated for the passage of an optical 
pulse into a semi-infinite sample. 


TAGS numbers: 42.50.Wk 


Theory of the radiation pressure on magneto-dielectric materials 2 

1. Introduction 

At the heart of the problem of radiation pressure is the famous Abraham-Minkowski 
dilemma concerning the correct form of the electromagnetic momentum in a material 
medium m El El E]; a problem which, despite of its longevity, continues to attract 
attention [3 El la El E]. The resolution of this dilemma lies is the identihcation of the 
two momenta, due to Abraham and Minkowski, with the kinetic and canonical momenta 
of the light, respectively [lO]. It serves to indicate, moreover, why the different momenta 
are apparent in different physical situations [9] . In this paper we shall be concerned with 
manifestations of optical momentum in radiation pressure on media and, in particular, 
on magneto-dielectric media. 

Most existing work on the theory of radiation pressure, as reviewed in [3 El 0 El E] > 
treated non-magnetic materials but there has been recent progress in the partial 
determination of the effects of material magnetisation. The rise of interest in 
metamaterials and in particular those with negative refractive index adds urgency to 
addressing this point m- 

We consider magneto-dielectric materials in which both the electric permittivity e 
and the magnetic permeability /i are isotropic functions of the angular frequency u. The 
quantum theory of the electromagnetic held for such media is now well-developed, with 
analyses in the literature based on Green’s functions US and also on Hopheld- 
like models based on coupling to a harmonic polarisation and magnetisation held 
[H113 [13 HU . The approach we shall adopt is to work with the elementary excitations 
within the medium, which are polarition coupled-modes of the electromagnetic held 
with the electric and magnetic resonances. We develop the associated classical linear- 
response theory and use this to determine the electric and magnetic held-huctnation 
spectra. The corresponding quantised held operators are expressed in terms of the 
polariton creation and destruction operators and we show that these satisfy the required 
standard commutation relations. We use these to construct the Minkowski and Abraham 
momentum operators and also the corresponding angular momenta, and show how the 
commutators of these with the vector potential facilitate the interpretation of the rival 
momenta [13 . The known form of the magnetic Lorentz force is conhrmed and used to 
explore the momentum transfer from light to a half-space sample and so provide the 
extension to permeable medium of earlier work on dielectrics |18] . 

2. Classical theory 

It suffices for our purposes to ignore the ehects of absorption and so work with a 
real permittivity and permeability. These arise naturally from the dispersion relations 
in a Hopheld-type model in which the electromagnetic held is coupled to harmonic 
polarization and magnetisation helds associated with the host medium HI. 
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The relative permittivity and permeability of a (lossless) magneto-dielectric at angular 
frequency u may be written in the simple forms 1201 


£(w) = n 


e=l 


/i(a;) = n 


^Te ~ 

^Im - 


m=l ^Tm - ’ 


( 1 ) 


where and are the numbers of electric and magnetic dipole resonances in 
the medium, associated with longitudinal and transverse frequencies indicated by the 
subscripts L and T. It follows that the values of the permittivity and permeability at 
zero and inhnite frequencies are 

ne 2 

<5(0) =nd^. <5(«5) = 1 

e=l ^Te 

TLm / 

A(o)=nT^. (2) 

m=l ^Tm 


in accord with the generalized Lyddane-Sachs-Teller relation [21]. The phase refractive 
index is 


Vp(^) = (3) 

as usual |22|. 

We consider a magneto-dielectric medium in the absence of free charges and 
currents, so that our electric and magnetic fields are governed by Maxwell’s equations 
in the form 


V-D 

V-B 

V X E 

V X H 


0 

0 

aB 

dt 

dB 

'W' 


(4) 


The electromagnetic helds described by these propagate as transverse plane-waves, with 
an evolution described by the complex factor exp(ik-r —icat), with k = |k| and u related 
by the appropriate dispersion relation. The four complex electric and magnetic helds 
are related by 122 ] 


D(a;) = £oE(a;) + P(a;) = 

B(w) = /io[H(a;) + M(a;)] = iaoia{u)B.{u ), (5) 


where P and M are, respectively, the medium’s polarization and magnetization. 
Henceforth, for the sake of brevity, we omit explicit reference to the frequency from 
our helds, permittivities and permeabilities. 
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Now suppose that external stimuli p and m with the same frequency are applied 
parallel to E and H respectively, so that the relations ([5]) become 

D = EqeE p 

B =/io/iH +/iom. (6) 

It is convenient to introduce 6-component held and stimulus variables dehned to be 
f = (E,ZoH) 

s = l/(p,m/c) (7) 

where Zq = \Jpio/^o is the usual free-space impedance and E is a suitably chosen sample 
volume. The energy of interaction between the held and the stimulus components is 
then 

6 

- EM = -^(E-P + hoH-m). 
j=i 

The solutions for the held components obtained by substitution of 
equations can be written 


( 8 ) 

into Maxwell’s 


f. = T.r 




(9) 


i=i 


where the linear response matrix is 

/ 


T = 


£oEDen 

V 

with the common denominator 

Den = , 



0 

0 

0 
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uckx 
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—ojckx 

0 

0 

—oockz 

OJCky 

U^E 

0 

0 

uckz 
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—uck^ 

0 

Oj'^E 

0 

—ojcky 

ojckx 

0 

0 

0 



■ ( 10 ) 


( 11 ) 


the zeroes of which give the dispersion relation for the held. The elements of the matrix 
T agree with and extend partial results obtained previously for non-magnetic media 


2.2. Field fluctuations 

The frequency and wave-vector huctuation spectra at zero temperature are obtained 
from the Nyquist formula [23] 

(/i/Ek = -Im(T,,) , (12) 

71 

where the required imaginary parts occur automatically in magneto-dielectric models 
that include damping mechanisms or, in the limit of zero damping, by the inclusion of 
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a vanishingly small imaginary part in w. In the latter case, which applies to our model, 
we hnd 


The total fluctuations are obtained by integration 



Im 


Den 


-)■ 


71 


2kc^ 




POO POO ^ P 

ififi) = du{fifi)^ = J dHfifdco, 


(14) 


where the frequency u and the three-diemnsional (3D) wave vector k are taken as 
continuous and independent variables. 

Consider plane-wave propagation parallel to the z-axis in a sample of length L and 
cross-sectional area A. The frequency fluctuation spectra are then obtained from the 
one-dimensional (ID) version of flTTl) as 

fj Tj roc 

{fifi)uj = ^ dklm{Tij) . (15) 

The contributions of the component x- and ^/-polarised waves give 



with the use of (ITOll and (|T^ . The E and H spectra, of course, have the usual relative 
magnitude for a magneto-dielectric material. The total fluctuations are obtained by 
integration of the spectra over u as in flTTll . 


2.3. Polariton modes 

The combined excitation modes of the material dipoles and the electromagnetic held 
are the polaritons [25]. The transverse polaritons for a material of cubic symmetry, as 
assumed here, are twofold degenerate corresponding to the two independent polarisation 
directions. Their dispersion relation is determined by the poles in the linear response 
function, the components of which are given in ffTOj) . The vanishing denominator (Den) 
gives 

, (17) 

which is the desired transverse dispersion relation. For our relative permittivity and 
permeability, there are rie + Um + 1 transverse polariton frequencies for each wave vector 
and these are independent of the direction of k. They jointly involve all of the resonant 
frequencies in the electric permittivity and the magnetic permeability. It is convenient 
to enumerate the transverse polariton branches by a discrete index u, and there is 
no overlap in frequency oj^u between the different branches. There are also 
longitudinal modes at frequencies cjLe and independent of the wave vector, but 
these will concern us no further. 
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Subsequent calculations involve the polariton group index rjg, defined in terms of 
the phase index rjp: 


dk 


(18) 


for Pp = ck/uj. The two refractive indices satisfy a satisfy a variety of sum rules over 
the polariton branches, given by [26l |27l [28] 

i:f—) =1 (19) 

U Vw/ku 


and by 


E 

u 

E 


£ 


u \hphg/k« 


= 1 

= 1 , 


( 20 ) 


where all of the optical variables are evaluated at frequency a;k«, so the sums run over 
every frequency Wku that is a solution of the dispersion relation flT7|) for a given wave 
vector k. 

3. Quantum theory 

3.1. Field quantisation 

The polaritons are bosonic modes, quantised by standard methods in terms of creation 
and destruction operators with the commutation relation 

«kii, d\.'u\ = , (21) 

The Dirac and Kronecker delta functions have their usual properties. The 
electromagnetic field quantisation derived previously for dispersive dielectric media [30] 
can be adapted for our magneto-dielectric medium by insertion of polariton branch 
labels and by conversion to SI units and continuous wave vectors. It is convenient to 
write all of the field operators in the form 

A(r,f) = A^+^(r, f) + , (22) 

where the second term is the Hermitian conjugate of the first. The deduced form of the 
transverse (or Coulomb gauge) vector potential operator is then given by 


A(+)(r,t) = 


h 


1/2 


167r%o 


1/2 


E / f ) hk. (t) e'*^"-ek 


\uJVpVgJ 


(23) 


where 


akjt) ^ iiku exp(-ic.Jku^) ■ 


(24) 
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The unit polarisation vector is assumed to be the same for all of the polariton branches 
at wave vector k. The two transverse polarisations for each wave vector are not shown 
explicitly here but they are important for the simplihcation of mode summations lai 
and they need to be kept in mind. For non-magnetic materials, with /i = 1, the vector 
potential fl2^ agrees with equation (5) from [32] and also with equation (A12) in [33] 
when the summation over polariton branches is removed. 

The electric held, E, and magnetic hux density or induction, B, held operators are 
readily obtained from the vector potential as 


E<+>(r,() 



(25) 


and 


B<+>(r.«) = VxA<+>(r,«) 

= I 


f dki ^ ) ak„(t) 


AU.r 


V16A„.' ^ ' -M ...... . — kxe,. (26) 

The remaining held operators, the displacement held D and the magnetic held H, follow 
from the quantum equivalents of ([5]) as 


D^^^(r, t) = i 


1/2 


IbTT^en 


E/^k 


f ^^Vp \ 

\ % J 


1/2 


dku{t)e^^'''ek 


(27) 


kix 


and 


/ t \ 1/2 / .. \ 1/2 

H(+^(r,t) = / ( ) E / ^k (-) ak„(t)e*‘"“'k X ek . (28) 

Vl67r%oy 

The above four operator expressions are consistent with previous work [HD]- The held 
quantisation for a magneto-dielectric medium has also been performed in a quantum- 
mechanical linear-response approach [I211I3] that assumes arbitrary complex forms for 
the permittivity and permeability. The resulting expressions for the E and B operators 
contain the polariton denominator fllll) in their transverse parts. A theory along these 
lines has been further developed by a more microscopic approach in which the explicit 
forms of the electric and magnetic susceptibilities are derived [Hcg. 


3.2. Field commutation relations 


The validity of the held operators derived here (or at least their self-consistency) 
is demonstrated by the conhrmation that they satisfy the required equal-time 
commutation relations [311 |3l]. Thus it follows from forms of our operators, given 
above, that 




ih 

Ibvr^ 

ih 

Svr^ 








kw 
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( 29 ) 


with recognition that the two exponents in the first step give the same contribntions 
and a crncial use of the hrst sum rule from equation flT^ . The transverse delta 
function, 6:^{r — r'), in the hnal step [311 IMl |35] relies on the implicit summation of the 
contributions of the two transverse polarisations. The common time t is omitted from 
the held operators in the commutators here and subsequently as it does not appear in 
the hnal results. An alternative canonical commutator is that for the vector potential 
and the electric held: 


Aj(r), -eoEj{v') 


ih 


E/^k 


= - r') , 



gik.(r-r') 




kti 


(30) 


with the use of a sum rule from equation fl20|) . It follows from the operator form of the 
hrst relation in equation ([5]) that the vectors potential and the polarisation commute: 


E(r),Pj(r ) 


= 0 . 


(31) 


This is most satisfactory as the two operators are associated with properties of diherent 
physical entities: the electromagnetic held and the medium respectively. We note 
that the commutator in equation (12^ has been given previously [lO] for non-magnetic 
dielectric media. 

The remaining non-vanishing held commutators are 


4(r),4(r') 


Svr^eo 


E/^k 


/i 






kti 


h 


= i—€ijhVf,5{ir - r') 
^0 




h 


Stt^IIoEo 


I dk[^] 


.VpVg, 


kli 


h 


= i - eijhVf^6{r - r') 


A(r),B,(r') 


h 


Vp_\ ir(r-r') 




Stt^ 

= iheijhVf,6{r - r') 


E-ijhkh 


A(r),4(r') 


h 


87r^/i( 


E / ^k 


.VpVg, 


^ir-{r—r') 


^ijhkh 


ku 


= i—€ijhVh5{r-v'), (32) 

To 

where we have used the sum rules given in equations flT^ and (120]) . Here Cijh is the 
familiar permutation symbol [36l [37] and the repeated index h is summed over the three 
cartesian coordinates x, y and z. The hrst of these results generalises a commutator 
previously derived for helds in vacuo [3l]. Note that together these commutation 
relations require that the polarisation operator commutes with the magnetic held 
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operator and that the magnetisation also commntes with the electric held operator: 


=0 

■M,(r),4-(r')l =0, 




which is a physical conseqnence of the fact that the two pairs of operators in each 
commntator are associated with properties of different systems: the medinm and the 
electromagnetic held. 


3.3. Parallel beams: 3D to ID conversion 


For parallel light beams, it is sometimes convenient to work with held operators dehned 
for dependence on a single spatial coordinate z with a one-dimensional wave vector k. 
Thns, for a beam of cross-sectional area A, conversion from 3D to ID is achieved by 
making the snbstitntions 


dk 




47r^ 


A 

A 


dk 


^(k - k') ^ -^S{k - k') 






27r 

The vector potential operator 
this way to 




(34) 


for polarisation parallel to the x-axis is converted in 




n 


1/2 


47r£o^, 


fdk 


1/2 




ikz . 


(35) 


where x is the nnit vector in the x-direction. The orthogonal degenerate polariton 
modes give a vector potential parallel to y. The fonr held operators retain their forms 
given in eqnations fl25|l to fl28l) except for appropriate changes in the hrst sqnare-root 
factors and vector directions. The electric and displacement helds are parallel to the x- 
axis, while the magnetic held and indnction are parallel to the //-axis. The commntation 
relation for the creation and destrnction operators retains the form given in eqnation 
bnt with k replaced by k. For /i = 1 and a single-resonance dielectric, the vector 
potential fl3^ agrees with eqnation (3.25) in [28] and eqnation (9) in [26] . 

The single-coordinate vector potential can also be expressed as an integral over 
freqnency by means of the conversions 


J dk —y 

6{k - k') 


(Iku 



— 6{u — u') 
Vg 



(36) 
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There are no overlaps in frequency between the different twofold-degenerate polariton 
branches and the summation over u is accordingly removed from the vector potential, 
which becomes 

A^+\z,t) = J du d^{t)exp{iu7]pz/c) (37) 

where 


a^(t) - a^e . 


(38) 


This destruction operator and the associated creation operator satisfy the continuum 
commutation relation 




5{u} — u'). 


(39) 


The held operators obtained by conversion in this way of fl2^ to (1281) are 

1 /2 

E^~''\z,t) =i J exp (zwr^pz/c) X 

r / hi ] \ 

B^^\z,t) =i J exp (ia;77p^/c) y 

t)^+\z,t) = i J du (soef^^iiaoiay^'^d^it) exp {iuT]pz/c) x 

ii(+)[z,t) =i J du du:{t)exp{iuripz/c)y . (40) 

There is again agreement with previously dehned expressions [2811^ for /i = 1. We can 
use these held operators to calculate the vacuum huctuations and hnd 

(0|Ej(^,i)|0) = (0|E)W(^,i)B)l->(^,()|0> = /(^)‘''' 

= {0|Bj(z,i)|0) 

(0|Hj(^,i)|0) = {0\E)<p(z,t)Hyp(z,m = jdojffj 

= {0|H=(j,()|0), (41) 


in exact agreement with results obtained in section 12.21 

We shall make use of the ID held operators derived here to calculate the force 
exerted by a photon on a magneto-dielectric medium, but hrst return to the full 3D 
description to investigate the electromagnetic momentum. 


4. Momentum operators 
f.l. Minkowski and Abraham 

The much debated Abraham-Minkowski dilemma is most simply stated as a question: 
which of two eminently plausible momentum densities, D x B and c“^E x H, is the 
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true or preferred value mig? We amplify upon the answer given in HDl, making special 
reference to the effects of a magneto-dielectric medium. 

Two rival forms for the electromagnetic energy-momentum tensors were derived 
by Minkowski [1] and Abraham [21 El i. Their original formulations considered 
electromagnetic helds in moving bodies, but it suffices for our our purposes to set the 
material velocity equal to zero. These results continue to hold for the magneto-dielectric 
media of interest to us. The two formulations differ principally in their expressions for 
the electromagnetic momentum and we consider here the respective quantised versions 
of these. 

The Minkowski momentum is quite difficult to hnd in the paper [1], but it can be 
deduced trom his expressions for other quantities. Its quantised form, is represented by 
the operator 


G^{t) = 


u,u' 


(irD(r, t) x B(r, t) 
1/2 


'dk(^] 

V % /ku 


/i 


1/2 




ak«(t)aL W + 


(42) 


where equations 
operator is 


and (ITT)) have been used. The diagonal part of this momentum 


'^diag ~ 2 



T hkuhj^^j h\^ 


The fact that k = rj^uj/c leads us to identify a single-photon momentum 

hu 

PiM) = - 



(43) 


(44) 


with the individual polariton mode hu. This form of the Minkowski single-photon 
momentum has been derived previously [3H] for a single-resonance non-magnetic 
material. There are good reasons, however as we shall see below, not to assign this 
value to the Minkowski momentum. 

The corresponding form of the Abraham momentum operator, proportional to the 
Poynting vector for energy flow, is 

1 


G"(t) = 


~2h 


(irE(r, t) x H(r, t) 
1/2 


dh 


Ufa 


) {—) 


1/2 


ku' 


ak«(t)aL'W + 4«Wak«'(t) 


t 


+hkM(^)h_kii'(t) + 4«W«-k«'W]k, (45) 

which differs from the corresponding expression for G'^(f) only in the the two square- 
root factors that occur in the held operators fl25D and fl28D . The diagonal part of this 
momentum operator is 


= - V 

^diag 9 Z-^ 


dk 


.VpVg 


) 




(46) 
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so that the associated single-photon momentum is 

huj f l\ 

Pa = — — , 

which is the usual Abraham value. 

The subscript M appears in brackets in equation 
form for the Minkowski momentum, given by 
huj 

Pm = —Vp 


(47) 


because there is an alternative 


(48) 


is observed in experiments sufficiently accurate to distinguish between the phase and 
group refractive indices, particularly the submerged mirror measurements of Jones and 
Leslie [39]. Note that the difference between the two candidate momenta. p(M and pm 
can be very large and even, the case of media with a negative refractive index, point in 
opposite directions [iQ]. The resolution of the apparent conflict between the two forms 
of Minkowski momentum is discussed in section 14.31 


4-2. Vector potential-momentum commutators 

It is instructive to evaluate the commutators of the two momentum operators with the 
vector potential. We can do this either by using the expressions for the helds in terms 
of the polariton creation and destruction operators or, more directly, by using the held 
commutation relations (129|1 and (IHOjl . together with the fact that the vector potential 
commutes with both the magnetic held and the induction. 

For the Minkowski momentum we hnd 

Ai{r),Gf = Cjki j dr' Ai{r),Dkir') Bfr') 

= -ihcjki J dr'(5,4(r - r')A(r') 

= -ihejkieimn J dr'6lj,{r - r')V^An{r') 

= - ihVjAir) +ihj dr'd^^ir - r ) 

= - ihVjAfr ), (49) 

where we have used the summation convention so that repeated indices are summed 
over the three cartesian directions. It should also be noted that the remaining helds 
will have the same form of commutation relation with the Minkowski momentum, for 
example for the electric held we have 

E,(r),Gf] =-thV,E,{r), (50) 

This follows directly from the relationships between these helds and the vector potential 
together with the fact that the vector potential commutes with the polarisation and the 
magnetisation. 
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For the Abraham momentum we can exploit our calculation for the Minkowski 
momentum to hnd the commutator 


Ai{r),Gf =ejki J dr' Ai{r), eoEk{r') fioHi{r') 

= Cjki j dr' Ai{r),eoEk{r') [Bfr') - 

= - ihVjAiir) + ihejki [ dr'S^f.{r - r')p,oMi{r') . 


(51) 


The second term, with its integration over the magnetisation, means that the 
commutator depends on both a held property, the vector potential, and a medium 
property, the magnetisation uni- 


4 . 3 . Interpretation 


The identihcation of the Minkowski and Abraham photon momenta respectively with 
the electromagnetic canonical and kinetic momenta has been proposed in the past 
[38111111121113], but rigorously proven only more recently [3 [10]. We need add here 
only a few brief remarks. 

The commutation relation satished by the Minkowski momentum operator 
resembles the familiar canonical commutator of the particle momentum operator. 




(52) 


where F(r) is an arbitrary vector function of position. Thus, analogously to its particle 
counterpart, the Minkowski momentum operator for the electromagnetic held generates 
a spatial translation, in this case of the vector potential and the electric and magnetic 
helds. The operator therefore indeed represents the canonical momentum of the held 
and it is the observed momentum in experiments that measure the displacement of a 
body embedded in a material host, as has been seen for a mirror immersed in a dielectric 
liquid [ 39 ], for the transfer of momentum to charge carriers in the photon drag ehect 
H and in the recoil of an atom in a host gas 05]. The simple spatial derivative that 
occurs on the right of equation fiT9|) shows that the measured single-photon momentum 
should have the Minkowski form in (l48il and not that given in equation (01]). 

The kinetic momentum of a material body is the simple product of its mass and 
velocity. The form of the Abraham single-photon momentum in equation (071) is verihed 
by thought experiments of the Einstein-box variety [IHIIIZ]. These use the principle of 
uniform motion of the centre of mass-energy as a single-photon pulse passes through 
a transparent dielectric slab and they reliably produce the Abraham momentum. The 
calculations remain valid with no essential modihcations when the slab is made from a 
magneto-dielectric material. 

More detailed analyses of the coupled material and electromagnetic momenta [3 HQ] 
show that the total momentum is unique but that this can be formed as the sum of 
alternative material and electromagnetic held contributions 


canonical 

medium 


= V 


kinetic 

medium 


+ G-‘ 


(53) 
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where the 'P operators represent the collective momenta of all the electric and magnetic 
dipoles that constitute the medium. The total momentum is the same for both 
the canonical and kinetic varieties, both being conserved in the interactions between 
electromagnetic waves and material media. 


4 . 4 . Angular momentum 


The electromagnetic held carries not only energy and linear momentum but also angular 
momentum and it is natural to introduce angular momenta derived from the Minkowski 
and Abraham momenta in the forms 


= / dr r X 


Dfr) X Bfr) 


= 


dr r X 


Efr) X Hfrl 


(54) 


A careful analysis of a light beam carrying angular momentum entering a dielectric 
medium shows that, in contrast with the linear momentum, the Minkowski angular 
momentum is the same inside and outside the medium, but that the Abraham angular 
momentum is reduced in comparison to its free-space value by the product of the phase 
and group indices |1H]. The analogue of the Einstein-box argument suggests that light 
carrying angular momentum entering a medium exerts a torque on it, inducing a rotation 
on propagation through it. An object imbedded in the host, however, may be expected 
to experience the influence of the same angular momentum as in free space and, indeed, 
this is what is seen in experiment |19]. 

The canonical or Minkowski angular momentum should be expected to induce a 
rotation of the electromagnetic helds, which requires both a rotation of the coordinate 
and also of the direction of the field. The requirement to provide both of these 
transformations provides a stringent test of the identification of the Minkowski and 
canonical momenta. It is convenient to hrst rewrite the Minkowski linear momentum 
density in a new form: 


D X B 


2 


DjViAj - DjVjAi 
DjWiAj - Vj [DjAi^ , 


(55) 


where we have used the hrst Maxwell equation, 'VjDj = 0. We can insert this form 
into our expression for and, on performing an integration by parts and discarding a 
physically-unimportant boundary term we hnd 


J^= dr 


Dj{r X 'V)Aj + D X A 


(56) 


which, in the absence of the medium reduces to the form obtained by Darwin [501 1^ - 
It is tempting, even natural, to associate the two contributions in the integrand with 
the orbital and spin angular momentum components of the total angular momentum. 
This is indeed reasonable, but it should be noted that neither part alone is a true angular 
momentum [52l [53l [M]- It is instructive to consider the commutation relation with a 
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single component of the angular momentum and so consider the operator 0 ■ J“: 

I ^-L 


A(r),0-J^' 

= — ih 

0 ■ (r X V)A 

^ + in{e 


= — ih 

0 ■ (r X V)A 

-ex A 


(67) 


The orbital and spin parts rotate, as far as is possible given the constraints of 
transversality, the amplitude and direction of the potential |52l|53llS]. The combination 
of both of these gives the required transformation. The commutator (l57|) gives the hrst 
order rotation of the vector potential about an axis parallel to 6 through the small angle 
6 , as the canonical angular momentum should. 


5. Magnetic Lorentz force 

It remains to determine the radiation pressure due to a light held on our magneto¬ 
dielectric medium. To complete this task we adopt the method used previously of 
evaluating the force exerted by a single-photon plane-wave pulse normally incident on 
the medium [m EH]. Before we can complete the calculation, however, we need to 
determine a suitable form for the electromagnetic force density. 

5.1. Heaviside-Larmor symmetry 

Maxwell’s equations in the absence of free charges and currents (|1|) exhibit the so-called 
Heaviside-Larmor symmetry [561 EZ], with the forms that are invariant under rotational 
duality transformations given by 122] 

E = E' cos f + ZqH' sin ^ 

H = H' cos sine 

D = D' cos e + sin ^ 

B = B' cose — sine , (58) 

for any value of f and where Zg is again the impedance of free space. It is readily 
verihed that the four Maxwell equations are converted to the same set of equations in 
the primed helds. The various physical properties of the electromagnetic held must also 
be unchanged by the transformations [HE]. We note, in particular, that the Minkowski 
and Abraham momentum operators, and given in equations and fH5|) and 
also the usual expressions for the electromagnetic held energy density and Poynting 
vector are all invariant under the transformation fl58l) . 

The standard form of the Lorentz force law in a non-magnetic dielectric is [59] 

fL = (p. v)E + /ioP X H, (59) 

with terms proportional to the electric polarisation and its time derivativ^. For a 
magneto-dielectric medium we need to add the force due to the magnetisation and to 

f This is usually written in terms of the magnetic induction as [59] 
fL = (p . V) E + P X B . 

In a magnetic medium, however, we need to distinguish between /ioH and B. That it should be the 
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do so in a manner that gives a force density that is invariant under the Heaviside- 
Larmor transformation. It follows from the transformation fl58|l that the polarisation 
and magnetisation are similarly transformed: 

P = P'cos.^ + c“^M'sin,^ 

M = M'cos ^ — cP'sin ^ . (60) 

The required form of the force density, satisfying the Heaviside-Larmor symmetry is 

fEL = (p . V) E + /ioP X H + /To (M ■ V) H - eoiaoM x E . (61) 

The invariance of this expression under the Heavisde-Larmor transformation is easily 
shown. This form of the force density was derived by Einstein and Laub over 100 
years ago [60] and there have since been several independent re-derivations of it 
[6T1 inn [63l mu EH ESJ- The £nal term in equation fl6T]) has been given special attention 
in IS7I, where it is treated as a manifestation of the so-called ‘hidden momentum’ given 
by £ohoM X E. This line of thought has attracted a series of publications, with several 
listed on page 616 of [22] • Omission of the hnal term leads to difhculties, not the 
least of which is the identification of a momentum density that does not satisfy the 
Heaviside-Larmor symmetry [68l l69j^ . The above derivation of the Einstein-Laub force 
density shows how the magnetic terms follow from the polarisation terms by simple 
symmetry arguments and so provides a new perspective on the complete force density. 
It is interesting to note, moreover, that the force density appropriate for a dielectric 
medium fl5^ may be obtained by consideration of the action on the individual dipoles 
making up the medium [59] and that the most direct way to arrive at the Einstein-Laub 
force density is to obtain the magnetisation part by treating a collection of Gilbertian 
magnetic dipoles m- 

It is also shown in the original paper [60] that the classical form of the Abraham 
momentum and the classical force density satisfy the conservation condition 

= - J dr f^^(r, t). (62) 

The integration is taken over all space and the relation is valid for helds that vanish 
at inhnity. This equality of the rate of change of the Abraham momentum of the light 
to minus the total Einstein-Laub force on the medium, or rate of change of material 
momentum, is as expected on physical grounds and further underlines the identihcation 
of the Abraham momentum with the kinetic momentum of the light mm- Equation 
fl6^ is simply an expression of the conservation of total momentum. 

former that appears in the force density follows on consideration of the screening effect of surrounding 
magnetic dipoles in the medium, in much the same way as electric dipoles screen the electric field in a 
medium. 

§ It is by no means straightforward to obtain the Einstein-Laub force density, in particular the final 
hidden-momentm related term, from the microscopic Lorentz force law and it has been suggested, for 
this reason, that the latter is incorrect [20]. A relativistic treatment, however, reveals that the required 
hidden-momentum contribution arises quite naturally from the Lorentz force law (ziiizaizsiiiiiza. 
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We consider a single-photon pnlse normally incident from free space at z < 0 on the 
flat snrface of a semi-inhnite magneto-delectric that hlls the half space z > 0. The 
pnlse is assnmed to have a narrow range of freqnencies centred on uq. Its amplitnde 
and reflection coefficients, the same as in classical theory, are m 

V^- 1 


R= - 


T = 


\J^/t + 1 
2 


(63) 


! s / p, + 1 

where all of the optical parameters depend on the freqnency. The damping parameters 
in the imaginary parts of £ and /r, as they occnr in R and T, are assnmed negligible bnt 
they shonld be snfficient for the attennation length 

( 64 ) 

2 a; Im(y^) 

to give complete absorption of the pnlse over its semi-inhnite propagation distance in 
the medinm. 

The momentnm transfer to the medinm as a whole is entirely determined by the 
free-space single-photon momenta, before and after rehection of the pnlse, as 


(l + R^) 


huo 



e + p huo 
2t]p c 


( 66 ) 


with the small imaginary parts of the optical parameters again ignored. It is often 
nsefnl to re-express the momentnm transfer in terms of a single transmitted photon 
with energy huo at z = 0~^. This qnantity is obtained by removal of the transmitted 
fraction of the pnlse energy, given by the bracketed term on the right of fl65|l . as 

£ -\~ p hojQ 1 


Ptotal 


2% 



( 66 ) 


in agreement with previons work [611 [77]. The remainder of the section considers the 
separation of this total transfer of momentnm to the magneto-dielectric into its snrface 
and bnlk contribntions. 

Previons calculations [iHl [55] of the radiation pressure on a semi-inhnite dielectric 
were made by evaluation of the Lorentz force on the material. This method is generalised 
here for a magneto-dielectric medium. For a transverse plane-wave pulse propagated 
parallel to the ^(-axis, with electric and magnetic helds parallel to the x and y axes 
respectively, the relevant component of the Einstein-Laub force from fl6T|) has the 
quantised form 

1 dM 


fEL _ TT 


dt 


-E 




( 67 ) 
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The ID field operators from section given in fl4Up . are appropriate here. The single¬ 
photon pulse is represented by the state 

|1) = y da;^(a;)a'^(a;)|0), (68) 

where |0) is the vacuum state. This single-photon state is normalised if we require our 
function to satisfy 

J duj= 1. (69) 

The states, from their construction, satisfy 




(70) 


A simple choice for the pulse amplitude is 


with c/L uo- The narrowness of the spectrum of this pulse means that u can often 
be set equal to Uq. 

The radiation pressure of the magneto-dielectric is obtained by evaluation of the 
expectation value of the Einstein-Laub force (|67|l for the single-photon pulse. The 
normal-order part of the force operator, indicated by colons, is used to eliminate 
unwanted vacuum contributions. A calculation similar to that carried out in [I8] for a 
dielectric medium gives 





/ y 



1 27rc2 y 

exp 

4c2 


{l\ : fj\z,t) \1) 



(72) 


where the real e, /i, and their derivatives in the group velocity are evaluated at frequency 
Uq and their relatively small imaginary parts survive only in the attenuation length i. 
It is readily verified by integration over t then z that this expression regenerates the 
total momentum transfer in equation fl66p . The force on the entire material at time t is 


{l\-.Ff\t)-.\l) = Af dz{l\-.ff^{zA)-.\l) 


EL/ 


(- 


huo 


X 


2 


^A/-erfc 


V 

/- ( ct L 

x/2 -- + 


%(£ + //) -277p 


ApU I 

2cH‘^ 


exp 


ct 


exp 


L2 


S^ip) 


, (73) 


with appropriate approximations neglecting small terms in the exponents for the long 
attenuation-length regime with L i. This expression, with error function and 
exponential contributions, has the same overall structure as found in other radiation 
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pressure problems [TKl 155] . The two terms in the large bracket are respectively the bulk 
and surface contributions, with time-integrated values 



(74) 


bulk 


surface 


This is again in agreement with equation (IM]) and it also reduces to equation (5.21) of 
[T8] for a non-magnetic material, where the momentum transfer can be written entirely 
in terms of the phase and group indices. The simple Abraham photon momentum again 
represents that available for transfer to the bulk material, once the transmitted part 
of the pulse has cleared the surface, while the more complicated surface momentum 
transfer depends on both £ and /i, together with their functional forms embodied in the 
phase and group refractive indices. 

We conclude this section by noting that the Heaviside-Larmor symmetry retains a 
presence in all of the forces and force densities obtained here, in that their forms are 
unchanged if we interchange, everywhere, the relative permittivity and permeability. 

6. Conclusion 

Much of the content of the paper presents the generalisations to magneto-dielectrics 
of results previously established for non-magnetic materials with /i = 1. We believe 
that the classical linear response theory in section 2 is novel; it allows, in particular, 
direct calculation of the electric and magnetic held-fluctuation spectra. The elementary 
excitations for a medium with multiple electric and magnetic resonances are the 
polaritons, whose phase and group velocities obey generalised sum rules for magneto¬ 
dielectrics [29] . 

The quantum theory in section 3 introduces electromagnetic held operators based on 
the multiple-branch polariton creation and destruction operators. It is shown that the 
vacuum huctuations of the quantised electric and magnetic helds reproduce the spectra 
obtained from our classical linear-response theory. The generalised held operators are 
shown to satisfy the same required canonical commutation relations as their simpler 
counterparts that hold in vacuum. 

The Minkowski and Abraham electromagnetic momentum operators are introduced 
in section 4 and their associated single-photon momenta are identihed. The 
commutators of these momentum and the vector-potential operators, previously 
calculated, rely on the canonical commutation relations and, through these, on the 
polariton sum rules. An extension to angular momentum, both canonical and kinetic, 
is achieved by introducing angular momentum densities that are the cross product of 
the position and the Minkowski and Abraham momentum densities respectively. 

Throughout our work we are guided by the Heaviside-Larmor symmetry between 
electric and magnetic helds. We show, in section 5, that application of this symmetry 
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leads directly to the Einstein-Laub force density [BU]. Our final result identifies the 
surface and bulk contributions in the force on a semi-inhnite magneto-dielectric for the 
transmission of a single-photon pulse through its surface. 
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